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Current spacecraft missions such as Wind and ACE can be used to determine magnetic correlation 
functions in the solar wind. Data sets from these missions can, in principle, also be used to compute 
so-called Eulerian correlation functions. These temporal correlations are essential for understanding 
the dynamics of solar wind turbulence. In the current article we calculate these dynamical cor- 
relations by using well-established methods. These results are very useful for a comparison with 
Eulerian correlations obtained from space craft missions. 
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O ■ I. INTRODUCTION 

(N ; 

In the past idealized models for turbulent fluctuations which can be found in the solar wind plasma or in the 
• interstellar medium have been proposed (e.g. Matthaeus et al. 1995). We are concerned with statistically axisymmetric 
models of magnetostatic fluctuations 5B{x) that are transverse to a uniform mean magnetic field Bq. If solar wind 
turbulence is considered, the mean field might be identified with the magnetic field of the Sun. The total magnetic 
\ field is a superposition of this mean field and the fluctuations B(x) = i?o + SB{x). Whereas we usually approximate 
^ ■ the mean field by a constant field aligned parallel to the z— axis {Bq = BqCz), the turbulent contribution has to be 
, replaced by turbulence models. Some prominant examples are slab, 2D, and two component models that include both 
' slab and 2D contributions (e.g. Matthaeus et al. 1990). 

There are recent spacecraft measurements of magnetic correlations in the solar wind (see e.g. Matthaeus et al. 2005, 
Dasso et al. 2007). Such measurements are very interesting and important since they allow an improved understanding 
■ of turbulence. For instance, characteristic length scales of turbulence such as the correlation length, the bendover 
scale, and the dissipation scale can be obtained from such observations. Also the investigation of spectral anisotropy 
by using data from different space craft missions such as Wind and ACE is possible. These properties of solar wind 
'" p^ ^ turbulence are very important for several investigations (heating and damping of the solar wind plasma, transport of 
I charged cosmic rays). A further important turbulence property is the turbulence dynamics (the time dependence of 
O I the stochastic magnetic fields). In principle, data sets from Wind and ACE can also be used to compute dynamical 
^ I correlation functions to explore the turbulence dynamics. 

^ ' In a recent article (Shalchi 2008) magnetic correlation functions were computed analytically. Such analytical forms 
• • of magnetic correlations complement data analysis results such as Matthaeus et al. (2005) and Dasso et al. (2007). 
Since we expect that future data analysis work will also allow the investigation of temporal correlation functions, we 
explore theoretically (numerically and analytically) the forms of these Eulerian correlations. These results can be 
5J] ' compared with data analysis results as soon as they are available. 
Ci The organization of the paper is as follows: in section 2 we define and discuss the basic parameters which are useful 

for describing turbulence. Furthermore, we explain the slab, the 2D, and the slab/2D compositel model. In section 
3 we review different models for the turbulence dynamics. In section 4 we compute Eulerian correlation functions 
numerically and analytically. In section 5 the results of this article are summarized. 
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II. GENERAL REMARKS - SETTING 



A. The turbulence correlation function 



The key function in turbulence theory is the two-point-two-time correlation tensor. For homogenous turbulence its 
components are 

Rimix,t) = (^dBiix,t)SB*^{6,0)) . (1) 
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The brackets (...) used here denote the ensemble average. It is convenient to introduce the correlation tensor in the 
fc— space. By using the Fourier representation 



SBi{x,t) = J (fk 5Bi{k,t)e'^-^ 

we find 

Rim{x,t) = j d'k j d'k i5Bi{k,t)5B*^{k ,0))e''^\ 
For homogenous turbulence we have 

'5Bi{k, t)6B*„^{k , 0) ) = Pi,^{k, t)5{k - k ) 



(2) 



(3) 



(4) 



with the correlation tensor in the /c— space Pim{k^t). By assuming the same temporal behaviour of all tensor compo- 
nents, we have 



Phn{k,t) = Pi„,{k)V{k,t) 

with the dynamical correlation funtion r(fc,t). Eq. Q becomes than 

Rim{x, t) - / d'k PUk)T{k, i)e^^-(^) 



(5) 



(6) 



with the magnetostatic tensor Pim{k) = (6Bi{k)6B*^{k)). 



B. The two-component turbulence model 

In this paragraph we discuss the static tensor Pim{k) defined in Eq. ([6]). Matthaeus & Smith (1981) have investigated 
axisymmetric turbulence and derived a general form of Pimik) for this special case. In our case the symmetry-axis 
has to be identified with the axis of the uniform mean magnetic field Bq = BqCz- For most applications (e.g. plasma 
containment devices, interplanetary medium) the condition of axisymmetry should be well satisfied. Furthermore, we 
neglect magnetic helicity and we assume that the parallel component of the turbulent fields is zero or negligible small 
{6Bz = 0). In this case the correlation tensor has the form 



Plrnik) ^ A{k\\,ki_) 



Sim — 



kl fcrjj 



(7) 



and Piz — Pzm — 0. The function A(A:||, k±) is controlled by two turbulence properties: the turbulence geometry and 
the turbulence wave spectrum. The geometry describes how A(k\\ , k±) depends on the direction of the wave vector k 
with respect to the mean field. There are at least three established models for the turbulence geometry: 



1. The slab model: here we assume the form 



S{kj 



A'""'{k\i,ki_) ^ g'^'^'iki^) 

In this model the wave vectors are aligned parallel to the mean field (fc || Bq). 
2. The 2D model: here we replace A(fc||, k±) by 

,-5(^11) 



(8) 



A2^(fc||,fc^)=g^^(fc^) 



2D/ 



(9) 



In this model the wave vectors are aligned perpendicular to the mean field (fc _L Bq) and are therefore in a 
two-dimensional (2D) plane. 



3 



3. The slab/2D composite (or two-component) model: In reality the turbulent fields can depend on all three 
coordinates of space. A quasi-three-dimensional model is the so-called slab/2D composite model, where we 
assume a superposition of slab and 2D fluctuations: 6B!^°^^{x,y, z) = 5Bf°-^{z) + SBf^{x,y). Because of 
< 5Bf°-^{z)SB*'^^{x,y) >= 0, the correlation tensor has the form 

PLT'ik) = PiTik) + P?^ik). (10) 

In the composite model the total strength of the fluctuations is SB^ — SB^^^j^ + ^B2D- The composite model 
is often used to model solar wind turbulence. It was demonstrated by several authors (e.g. Bieber et al. 1994, 
1996) that 20% slab / 80% 2D should be realistic in the solar wind at 1 AU heliocentric distance. 

The wave spectrum describes the wave number dependence of A{k\\ , k\^. In the slab model the spectrum is described 
by the function g''^°'^{k\\) and in the 2D model by g'^^{k±). 

As demonstrated in Shalchi (2008), the combined correlation functions (defined as R± = Rxx + Ryy) for pure slab 
turbulence is given by 

/■OO 

Rf^^z) ^Stt dkn g^'^^k^i) cos{ki^z) (11) 
Jo 

and the correlation function for pure 2D is 

Rf{p) = 2n dk^g^''(kj_)Jo{k^p). (12) 

JQ 

Here z is the distance parallel with respect to the mean magnetic field and p denotes the distance in the perpendicular 
direction. To evaluate these formulas we have to specify the two wave spectra g'*'"^(fc||) and g^^{k±). 

C. The turbulence spectrum 

In a cosmic ray propagation study, Bieber et al. (1994) proposed spectra of the form 

^2D(fc^) ^ ^l^^SB-Ul + kiq^)--^ (13) 
TT 

with the inertial range spectral index 2i/, the two bendover length scales I slab and I2D, and the strength of the slab 
and the 2D fluctuations SB^i^^^ and SB^jj. By requiring normalization of the spectra 

SB^ = 5Bl + 5Bl + 5Bl = j d?k [p,,(fc) + Pyy{k) + P,,{k)\ (14) 

we find 

27^1X^7^- ^''^ 

By combining these spectra with Eqs. fTT]) and (fT2)l the slab correlation function 

Ri^'iz) = ^C{v)5Bl,,lsiab / dfc|i (1 + klll,,)-^ coB{k\\z) (16) 

Jo 

as well as the 2D correlation function 

/"OO 

i?f (p) = AC{u)5Blj,hD / (1 + klllj,)--Uk^p) (17) 

Jo 

can be calculated. In Eq. ([T7]) we have used the Bessel function Jo (a:). In Shalchi (2008) such calculations valid for 
magnetostatic turbulence are presented. 
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Model 






Magnetostatic model 

Damping model of dynamical turbulence 

Random sweeping model 

Undampled shear Alfven waves 

Undampled fast mode waves 

NADT model 


1 

^ — CtV A fell t 

COs(±UAfc||t) 

cos{vAk\\t) 
cos(±i;Afc||i)e~^='°''* 


1 

e ^ ^ 

^-(aVAk±t)^ 
1 

cos{vAk±t) 

g — f2nt 



TABLE I: Different models for the dynamical correlation function r{k,t). Here, va is the Alfven speed and a ts a parameter 
that allows to adjust the strength of dynamical effects. The parameters jsiat and of the NADT model are defined in Eq. 



D. Correlation functions for dynamical turbulence 

For dynamical turbulence the slab and the 2D correlation functions from Eqs. (fTTjl and (fT2|) become 

/>oo 

Rf^^z) = Stt / g"''^^(fcj|)cos(fc||z)r^'"''(/c||,t) 
Jo 

POO 

Rl^ip) = 2TT / dk^g^^{kMk±p)r'''{k^,t). 
Jo 

For the model spectrum defined in the previous paragraph these formulas become 

Rt\z,t) = 4C7(i^)5BL6^.;a6 / (l + fc2/2^^,)--cos(fc||z)r^''^''(fc||,i) 

= 4c{i^)SBiohD / dk^{i + kiqo)-'jo{k^p)r^''{k^,t). 



(18) 



(19) 



To evaluate these equations we have to specify the dynamical correlation functions r^''^''(fe|| , t) and T'^^{k±,t) which 
is done in the next section. 

III. DYNAMICAL TURBULENCE AND PLASMA WAVE PROPAGATION EFFECTS 

In the following, we discuss several models for the dynamical correlation function T{k,t). In Table |T1 different 
models for the dynamical correlation function are summarized and compared with each other. In the fc— space these 
models show a very different decorrelation in time. In the Section IV we compute the dynamical correlation functions 
in the configuration space for these different models 



A. Damping and random sweeping models 

One of the first authors discussing dynamical turbulence were Bieber et al. (1994). In their article, the authors 
proposed two models for the dynamical correlation function: 



r_DT(fc, i) = e (damping model of dynamical turbulence) 

TRs{kjt) — 6"^*/*'=*^'^^)^ (random sweeping model) 

with the correlation time scale tc{k). Bieber et al. (1994) estimated the correlation time as 

1 



tc{k) 



avA I k 



(20) 



(21) 



Here, va is the Alfven speed and a is a parameter which allows to adjust the strength of the dynamical effects, ranging 
from a = (magnetostatic turbulence) to a = 1 (strongly dynamical turbulence). Bieber et al. (1994) also suggested 
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that the parameter a could be interpreted as SB/Bq. In this case, the correlation time scale tc{k) becomes comparable 
to the eddy turnover time. Also, decorrelation effects related to plasma waves (see, e.g., Schlickeiser & Achatz 1993) 
can be achieved by expressing a through parameters such as the plasma (3 (for a definition see also Schlickeiser & 
Achatz 1993). The damping model was originally introduced for a particle scattering study in dynamcial turbulence 
by Bieber et al. (1994). In this model, the dynamical correlation function has an exponential form, whereas in the 
random sweeping model T(k, t) has a Gaussian form. 



B. Plasma wave turbulence 



Another prominent model is the plasma wave model which is discussed in Schlickeiser (2002). In this model, the 
dynamical correlation function has the form 

Tpw{k, t) = cos(wi)e-'^*. (22) 

Here, uj is the plasma wave dispersion relation, whereas 7 desribes plasma wave damping. Often, undamped plasma 
waves are considered, where T pw{k^t) = cos(wi) and the dynamical correlation function is a purely oszillating 
function. Prominent examples for different plasma waves are Shear Alfven waves, where lo = ±w^fc||, and fast 
magnetosonic waves, where m — VAk. 



C. The nonlinear anisotropic dynamical turbulence model 

Recently, an improved dynamical turbulence model, namely the nonlinear anisotropic dynamical turbulence 
(NADT) model, has been proposed by Shalchi et al. (2006). This model takes into account plasma wave propa- 
gation effects as well as dynamical turbulence effects. The NADT model was formulated for the slab/2D composite 
model, where, in general, we have the two different dynamical correlation functions r'*''^''(fc|| , t) and r^^{k±, t), namely 

r'"^(fc||,t) = cos{ujt)e-^''^' ' 

r2^(fcx,t) = e-^'°* (23) 

with 

2D n j 1 for k±_l2D < 1 

^ = ^ 1 {kj,^?/^ for kj,^ > 1. (24) 

and with the plasma wave dispersion relation of shear Alfven waves 

UJ = jvAk\\ i = ±l. (25) 

In Eq. the parameter f3 can be expressed by the strength of the 2D component 6B2d/Bq, the 2D bendover scale 
I2D, and the Alfven Speed va (see Shalchi et al. 2006): 

(26) 

ho i'o 

These forms of the temporal correlation function are discussed in more detail in Shalchi et al. (2006). They are based 
on the work of Shebahn (1983), Matthaeus et al. (1990), Tu & Marsch (1993), Oughton et al. (1994), Zhou et al. 
(2004), Oughton et al. (2006). In the current article we approximate 7^^ in Eq. ((24l) by 

^P{l + kj2Df' (27) 

for simplicity. 

The parameter j in Eq. (|25p is used to track the wave propagation direction (j = +1 is used for forward and j = —1 
for backward to the ambient magnetic field propagating Alfven waves) . A lot of studies have addressed the direction 
of propagation of Alfvenic turbulence, see, for instance, Bavassano (2003). In general, one would expect that, closer 
to the sun, most waves should propagate forward and, far away from the sun, the wave intensities should be equal 
for both directions. Most of the observations, which allow conclusions on space plasma and particle propagation 
properties, have been performed in the solar wind at 1 AU heliocentric distance. Thus, we can assume that all waves 
propagate forward, and we therefore set j = +1 in the current article. 
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IV. EULERIAN CORRELATION FUNCTIONS 



by 



To investigate the different models we calculate the (combined) single-point-two-time correlation function defined 



E^{t) := R^{x = 0,t) = {6Bi{0, t)SB*,{Q, 0)) 



(28) 



Since this function is of particular importance in understanding dynamical turbulence effects and interactions with 
energetic charged particles this function is also known as Eulerian correlation function Ej_ (t) . For the different models 
discussed in Section III, Eqs. (flQl) become 



and 





1 


magnetostatic model 




cos(Ta;) 


Alfven waves 


X < 




damping model 






random sweeping model 




cos(Ta;)e~'''^ 


NADT model 


= 4C{i^)6BId / 
./n 


dx (l + a;^)-'' 




1 


magnetostatic model 




1 


Alfven waves 


X < 


'slab 

— a , Tx 

e '2D 


damping model 




e '2D ' 


random sweeping model 






NADT model 



(29) 



(30) 



Here we have used the integral transformations x = k\^lsiah and x = k±l2D- Furthermore we used the dimensionless 
time 



and the parameter 



^ = V2 



V At 1 1 slab 
6B2D Llab 



Bq I2D 

In the following paragraphs we evaluate Eqs. (^5)1 and ([50]) numerically and analytically. 



(31) 
(32) 



A. Numerical calculation of Eulerian correlations 



In Figs. [Hand [2] the results for the slab correlation function and in Figs. [3] and [4] for the 2D correlation function 
are shown for the different dynamical turbulence models. To obtain these figures we have solved the integrals in Eqs. 
(|29p and pop numerically. For the damping model and the random sweeping model we used a = 1. Furthermore, 
we used I2D = O.llgiab as in previous articles based on the result of laboratory experiments such as Robinson & 
Rusbridge (1971). For the turbulence spectrum in the inertial range we employ a Kolmogorov (1941) behavior by 
setting 1/ = 5/6. 

As shown in Figs. [T]-|31 the Eulerian correlations obtained for the damping model and the random sweeping model 
are very similar. The results obtained by employing the NADT model are, however, quite different from the other 
models. 

B. Analytical calculation of Eulerian correlations 

Here we compute analytically the different Eulerian correlations. For magnetostatic (MS) turbulence we can use 

r°° 1 



0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 

t: = v. t/l , ^ 
A slab 



FIG. 1; The slab correlation function Rf"-^ {Q,t) / SBli^f, as a function of the time r = VAt/hiab- Shown are the results obtained 
for the Alfvenic plasma wave model (dotted line), the damping model of dynamical turbulence (dashed line), the random 
sweeping model (dash-dotted line), and the NADT model (solid line). 

and therefore 

which is the expected result. In the following paragraphs we investigate the different other turbulence models. 



1. Undampled shear Alfven waves 



In this case we can use (see, e.g., Shalchi 2008) 

1 /2\^/^"'' 

dx (1 + cos(rx) - ^ f - j i^i/2-.(T) (35) 



to derive 



Obviously the characteristic time scale for temporal decorrelation tc is 

^slab.Alf _ hlab ^^^-^ 
VA 

Following Shalchi (2008), the Modified Bcsscl function in Eqs. ([55]) and can be approximated for large arguments. 
We find for times much larger than the temporal decorrelation time 

^slab,Alf^^ » Q « (38) 

For the special case = 1 we obtain an exponential function. For the 2D Eulerian correlations we always have 
j^2D,Aif _ gj^2^^ since there are no wave propagation effects in the perpendicular direction for (undamped) 
Alfvenic plasma waves. 
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t: = v. t/l , , 
A slab 



FIG. 2; Same caption as in Fig. [T]but now as semi- log-plot. Clearly we can see that we find an exponential function for the 
Eulerian correlation if we employ the Alfven wave model (dotted line) or the NADT model (solid line). 



2. Damping model of dynamical turbulence 

For the damping model of dynamical turbulence (DT) the integrals in Eqs. ([29|) and ([30|) arc difficult to solve. The 
results can be found in the appendix. As shown there, the characteristic time scale for the temporal decorrelation is 



,DT 



for slab fluctuations 

avA 

^ for 2D fluctuations 

avA 



(39) 



For the case t ^ tc corresponding to OiT » 1 (i = slab, 2D) with 



O^slab 



0-20 



ho 



(40) 



we can easily compute the Eulerian correlation function approximatelly. For large a^r there is only a contribution to 
the exponential function in Eqs. (I29|) and (j30p for a; — ^ 0. Thus, we can approximate 



2 \—u—aiTX 



OiT 



(41) 



to obtain 



AC{v) j 5B^i^fJ,siab for slab fluctuations 



avAt 1 5B2£il2D for 2D fluctuations. 
For the damping model of dynamical turbulence the Eulerian correlation function tends to zero with - t-^ 



(42) 



3. Random sweeping model 



For the random sweeping model the analytical results can also be found in the appendix. As shown there we find 
the same temporal correlation time scale tc as for the damping model of dynamical turbulence (see Eq. p9p ). For 
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1.5 
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t: = v. t/l , ^ 
A slab 



3.5 



4.5 



FIG. 3; The 2D correlation function {0,t) /SB2D a^s a function of the time r = VAt/lsiab- Shown are the resuhs obtained 
for the Alfvenic plasma wave model (dotted line), the damping model of dynamical turbulence (dashed line), the random 
sweeping model (dash-dotted line), and the NADT model (solid line). The result for undampled Alfven waves corresponds to 
the magnetostatic model. 



time scale satisfying t we can use 



2a,; r 



to find 



^r,q, , \ SBl,,^^^^ for slab fluctuations 

^;r(t»t.) = 2V^C(.) £^ fo. 2D fluctuations. 



(43) 



(44) 



Obviously the results for the random sweeping model are very similar to the results obtained for the damping model. 
This conclusion based on analytical investigations agrees with the numerical results from Figs. [T]-[3) 



4. NADT model 



Here we have to distinguish between the slab and the 2D correlation function. For the slab function we can use the 
result for Alfven waves with an additional factor ea;p(— ^r). Therefore, wc find for late times 



r(^.-i/2) 



slab 



VAt 



'VAt(l + 0/hla.b 



(45) 



In this case there are two correlation time scales. The first is associated with the plasma wave (PW) propagation 
effects 



slab, NADT 
c,PW 



^slab 

VA 



and the second is associated with the dynamical turbulence (DT) effects. The latter correlation time is 



.NADT 



slab 



1 Bo I 



2D 



"'""^ VA^ V2SB2D VA 



(46) 



(47) 
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0.5 1 1.5 2 2.5 3 3.5 4 4.5 

t: = v. t/l , , 
A slab 



FIG. 4: Same caption as in Fig. |3]but now as semi- log-plot. Clearly we can see that we find an exponential function for the 
Eulerian correlation if we employ the NADT model (solid line). 

For the 2D fluctuations the situation is more complicated and, thus, the analytical calculations can be found in the 
appendix. As demonstrated there the correlation time scale is given by Eq. (|47p . The behavior of the Eulerian 
correlation function for late time {t 3> t^or'^) exponential function 

El^ ^AC{iy)SBlae~''^'^^^'""'. (48) 
This exponential result agrees with our numerical findings visualized in Figs. [3]and[4l 

V. SUMMARY AND CONCLUSION 

In this article we have calculated and discussed Eulerian correlation functions. The motivation for this work are 
recent articles from Matthaeus et al. (2005) and Dasso et al. (2007). In these papers it was demonstrated, that 
magnetic correlation functions can be obtain from spacecrafts measurements (ACE and Wind) . We expect that from 
such observations also Eulerian correlations can be obtained. In the current article we computed analytically and 
numerically these correlations. These theoretical results are very useful for a comparison with data obtained from 
ACE and Wind. 

We have employed several standard models for solar wind turbulence dynamics, namely the (undamped and 
Alfvenic) plasma wave model, the damping model of dynamical turbulence, the random sweeping model, and the 
nonlinear anisotropic dynamical turbulence (NADT) model. All these model are combined with a two-component 
model and a standard form of the turbulence wave spectrum. As shown, we find very similar Eulerian correlations 
for the damping model and the random sweeping model. Therefore, we expect that in a comparison between these 
models and spacecraft data, one can not decide which of these models is more realistic. The NADT model presented 
in Shalchi et al. (2006), however, provides different results in comparison to these previous models. In table [TTl we 
have compared the different correlation time scale derived in this article for the different models. 

By comparing the results of this article with spacecraft measurements, we can find out whether modern models like 
the NADT model is realistic or not. This would be very useful for testing our understanding of turbulence. Some 
results of this article, such as Eqs. are quite general and can easily be applied for other turbulence models (e.g. 
other wave spectra). 
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Model 


isLab 
'•c 


,'2U 
'•c 


Magnetostatic model 

Undampled shear Alfven waves 

Damping model of dynamical turbulence 

Random sweeping model 

NADT model (plasma wave effects) 

NADT model (dyn. turbulence effects) 


CXD 

hlab 

"A 
hlab 
avA 
Islab 
aVA 
Islab 

■"A 

1 So l2n 


CXD 
OO 
ho 

avA 

hp 

aVA 

no effect 

1 So ho 


•J2 SB2D 


•J2 SB2D ^A 



TABLE If: Comparison for the different correlation time scale found in this article. For the damping model of dynamical 
turbulence and the random sweeping model we found the same correlation times. For the NADT model there are two correlation 
times, one scale for the plasma wave propagation effects and one scale for the dynamical turbulence effects. 
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APPENDIX A: EXACT ANALYTICAL FORMS OF EULERIAN CORRELATION FUNCTIONS 

By using Abramowitz & Stegun (1974) and Gradshteyn & Ryzhik (2000) we can compute analytically the different 
Eulerian correlation functions defined in Eqs. (|29p and (I30p . For the plasma wave model the results are given in the 
main part of the paper. 

1. Damping model of dynamical turbulence 

For the damping model of dynamical turbulence we can use 

rfa; (l + x2)-'^e-»--- = ^^^11 

X [2Ji,_i/2{atT)Csc{2'Kv) - Ji/2-„{aiT)Sec{'Kiy) + Csc{Trv)Hi/2-^{aiT)] . (Al) 







Here we used Bessel functions Jn{z) and the Struve function iJ„(z). By making use of this result we find for the 
Eulerian correlation function 



2-l/2-.^3/2(^^^).-l/2 



X [2J^_i/2{aiT)C sc(2ny) - Ji/2-i,{aiT)Sec(TTv) + Csc(7riy)iJi/2_y (a^r)] (A2) 

with i = slab, 2D and 

dslab = a 

Islab / A o\ 

a2D = a- — (A3) 
ho 

for the damping model of dynamical turbulence. 



2. Random sweeping model 

For the random sweeping model we can employ 



' dx (1 + x2)-e-(--)^ ^^u(y--., ] (A4) 
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with the confluent hypergeometric function U{a,b,z). By employing this result we find for the Eulerian correlation 
function 

E^f'it) = 2y^C{,.)SB^U J - (A5) 

for the random sweeping model. 



3. NADT model 



For the NADT model we only have to explore the 2D fluctuations (the slab result is trivial and discussed in the 
main part of the this paper). Eq. (|30p can be rewritten as 



dx e-«^ 



dx x~'^''e-^^ 



The first integral in trivial, the second one can be expressed by an exponential integral function E^{z): 



El"" = AC{v)5Bli, 



(A6) 



(A7) 



This is the final result for the Eulerian correlation function of the 2D fluctuations. To evaluate this expression for 
late times (^r 3> 1) we can approximate the exponential integral function by using 



i?3.-l/2(CT»l)«^ 



(A8) 



Here we assumed that the main contribution to the integral comes from the smallest values of x, namely x w 1. By 
combining Eq. (|A8p with Eq. (|A7p we find approximatelly 



E^^ w AC{v)5Bljje-^^ 



(A9) 



corresponding to an exponential behavior of the Eulerian correlation function. For the correlation time scale we find 
Tc = 1/^. A further discussion of these results can be found in the main part of the text. 
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